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Abstract: Stable isotopes in tree-ring cellulose provide important data in ecological, archaeological, and 
paleoenvironmental researches, thereby, the demand for stable isotope analyses is increasing rapidly. 
Simultaneous measurement of cellulose 61°C and 6!8O values from tree rings would reduce the cost of 
isotopic commodities and improve the analytical efficiency compared with conventional separate 
measurement. In this study, we compared the 813C and 6!8O values of tree-ring a-cellulose from Tianshan 
spruce (Picea schrenkiana) in an arid site in the drainage basin of the Urumqi River in Xinjiang of 
northwestern China based on separate and simultaneous measurements, using the combustion method (at 
1050°C) and the high-temperature pyrolysis method (at 1350°C and 1400°C). We verified the results of 
simultaneous measurement using the outputs from separate measurement and found that both methods 
(separate and simultaneous) produced similar 6!3C values. The two-point calibrated method improved the 
results (range and variation) of 6!3C and 6!8O values. The mean values, standard deviations, and trends of 
the tree-ring LC obtained by the combustion method were similar to those by the pyrolysis method 
followed by two-point calibration. The simultaneously measured 6'8O from the pyrolysis method at 
1400°C had a nearly constant offset with that the pyrolysis method at 1350°C due to isotopic-dependence 
on the reaction temperature. However, they showed similar variations in the time series. The climate 
responses inferred from simultaneously and separately measured 61°C and 6'8O did not differ between the 
two methods. The tree-ring BC and 6!8O values were negatively correlated with standardized 
precipitation evapotranspiration index from May to August. In addition, the 6'8O was significantly 
correlated with temperature (positive), precipitation (negative), and relative humidity (negative) from May 
to August. The tree-ring ÒC and 6!8O values determined simultaneously through the high-temperature 
pyrolysis method could produce acceptable and reliable stable isotope series. The simultaneous isotopic 
measurement can greatly reduce the cost and time requirement compared with the separate isotopic 
measurement. These results are consistent with the previous studies at humid sites, suggesting that the 
simultaneous determination of VC and 6'8O in tree-ring a-cellulose can be used in wide regions. 
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1 Introduction 


Tree-ring stable isotopes provide important data for analyses in dendroclimatology, 
dendroecology, and paleoecology (Schleser et al., 1999; McCarroll and Loader, 2004; Babst et al., 
2014). The measurement and interpretation of stable isotopes in tree-ring cellulose are 
increasingly influential and are rapidly improving the development of dendroclimatology, 
dendroecology, and paleoecology (McCarroll and Loader, 2004; Dawson and Siegwolf, 2007; 
Werner et al., 2012). These studies have greatly improved our understanding of the past 
environmental variation (Treydte et al., 2007; Leavitt et al., 2010; Liu et al., 2017) and of the 
plant responses to the current environment or the changing climate (Saurer et al., 2014; Altieri et 
al., 2015; Frank et al., 2015). Recently, stable isotopes records in tree rings have been 
incorporated into isotope-enabled Earth system models to evaluate the large-scale mechanisms of 
variation in photosynthesis and stomatal conductance under a changing climate and environment 
(Saurer et al., 2014; Keel et al., 2016; Keller et al., 2017). However, the veracity of the 
implementation of these models needs to be confirmed by a closer examination of regional and 
long-term isotopic chronologies in tree rings. 

The stable carbon isotope (8'C) and stable oxygen isotope (5!8O) are traditionally measured 
separately. Unfortunately, high cost of sample preparation coupled with time-consuming 
analytical procedures are required to obtain a high level of replication and an adequate spatial 
coverage which are required for robust paleoclimatology studies and are difficult to achieve 
(Leavitt et al., 2010). Obtaining stable isotope measurements based on a long-term chronology is 
difficult due to the limited quantity of wood available in narrow tree rings from arid regions 
(Wang et al., 2013) and from archaeological studies (such as subfossil wood) (Ktusek and 
Pawetczyk, 2014). Furthermore, when trends in two or more stable isotopes are viewed together, 
this can provide far greater insights into plant physiology and environmental responses than the 
stable isotopes that are considered separately (Loader et al., 2015). For example, the relationship 
between 8!°C and 5!8O reflects a balance between plant carbon and water metabolism (Babst et 
al., 2014; Liu et al., 2014). Tree-ring 5'°O mainly reflects the 5!8O of the source water and 
environmental humidity, whereas tree-ring 5!°C reflects the variations in photosynthetic rate and 
stomatal conductance (Farquhar et al., 1982; Roden et al., 2000; Roden and Farquhar, 2012). 
Therefore, the combination of 5'°C and 5'8O could capture more environmental signals 
simultaneously than a single isotope, and could improve our understanding of the plant's physical 
response to climate change and of the terrestrial carbon cycle in relation to climate change 
(Scheidegger et al., 2000; Babst et al., 2014). Thus, it would be beneficial to develop rapid and 
high-efficiency methods for a reliable stable isotope analysis by combining the two stable 
isotopes (i.e., 54°C and 5!80). 

To accomplish this goal, some researchers investigated the potential for simultaneous 
measurement of 5!3C and 5!8O from low- or high-temperature pyrolysis of cellulose (KnGller et 
al., 2005; Young et al., 2011; Woodley et al., 2012). The time and cost of such an approach would 
decrease substantially, since multiple isotopic ratios could be analyzed in the same time, without 
an extra expense. Several comparative studies have reported encouraging results with this 
approach (Saurer et al., 1995; Kndller et al., 2005; Young et al., 2011; Woodley et al., 2012). 
However, these studies were conducted in humid areas (Kndller et al., 2005; Young et al., 2011; 
Woodley et al., 2012). Few researches have been conducted in arid areas, where the range of 5'8O 
value is relative small (Xu, 2014). 
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The purposes of the present study were to confirm the applicability of simultaneous 
measurement of calibrated 8'°C and 5!8O in tree-ring a-cellulose from Tianshan spruce (Picea 
schrenkiana) in climatic and ecological studies and to determine whether it is possible to compare 
stable isotope values obtained from this method with values from separate measurement in a 
previous research conducted in an arid site in Xinjiang of northwestern China. First, we 
determined the relationships between the reaction temperature and the simultaneously measured 
SLC and 5!8O values to provide a reference for the comparison of stable isotope estimates based 
on the pyrolysis method at different reaction temperatures. In this analysis, we confirmed the 
viability and stability of the pyrolysis method for the simultaneous measurement of 51°C and 8'8O 
from tree-ring a-cellulose. Next, using previous data from separate measurements of 51°C 
(combustion analysis) and 6!8O (pyrolysis at 1350°C) in tree-ring a-cellulose from P. schrenkiana, 
we compared the values of lC and 5'8O from separate measurement with those from 
simultaneous measurement by the high-temperature (at 1400°C) pyrolysis of tree-ring o-cellulose. 
Finally, we evaluated whether the 8'°C and 6!8O values measured simultaneously by the 
high-temperature pyrolysis of tree-ring a-cellulose could be used for paleoenvironmental and 
paleoecological researches in arid regions. 


2 Data and methods 


2.1 Materials and experimental design 


We obtained tree-ring samples from P. schrenkiana growing in an arid site in the drainage basin of 
the Urumqi River in Xinjiang Uygur Autonomous Region of northwestern China. The tree-ring 
cores were obtained from mature trees using a 12-mm-diameter increment borer at breast height 
(about 1.3 m above the ground) in September 2008. We totally collected 40 cores from 20 mature 
trees (two cores per tree) and processed these cores according to procedures described by Stokes 
and Smiley (1968). The tree-ring width was measured and cross-dated using a Lintab device 
(LINTAB 6; Rinntech, Heidelberg, Germany) connecting with the Time Series Analysis and 
Presentation dendrochronological software. We selected 10 of these cores at a rate of one core per 
tree and pooled the wood materials for a given year to extract a-cellulose based on the method of 
Loader et al. (1997). The a-cellulose was extracted and then homogenized using an ultrasound 
machine (JY92-2D, Ningbo Scientz Biotechnology Company, Ningbo City, China) (Laumer et al., 
2009). The composite samples were first used for the separate analyses in our previous study (Xu 
et al., 2014) and the remaining samples were used for the simultaneous analyses in this study. 

The 5!°C was measured in 2010 using the conventional method of combustion at 1050°C in an 
elemental analyzer (Flash EA 1112, Thermo Fisher Scientific, Milan, Italy) coupled with a 
continuous-flow isotope ratio mass spectrometer (Delta Plus, Thermo Electron Corporation, 
Bremen, Germany) at the Key Laboratory of Western China's Environmental System, Ministry of 
Education, Lanzhou University, China. In the combustion analysis, the tin capsules packaged 
0.18—0.20 mg a-cellulose were transformed into CO2 in the presence of excess oxygen (Young et 
al., 2011). The analytical uncertainty was 0.05%o (=1o) for 6'3C. For 8'8O measurement, we 
loaded 0.14—0.16 mg of a-cellulose into silver capsules and determined the isotope ratio using a 
high-temperature conversion elemental analyzer (TC/EA; Thermo Electron Corporation, Bremen, 
Germany) at 1350°C coupled to a Finnigan MAT-253 mass spectrometer (Thermo Electron 
Corporation, Bremen, Germany) at the State Key Laboratory of Cryospheric Science, Northwest 
Institute of Eco-Environment and Resources, Chinese Academy of Sciences (CAS), Lanzhou, 
China. Hereafter, we refer to this 8'8O chronology as the "pyrolysis-1350" series. More details of 
the sample preparation and stable isotope measurements were provided in Xu et al. (2014). 

The simultaneous measurement of 5'°C and 8'8O was conducted using a TC/EA coupled to a 
Finnigan MAT-253 mass spectrometer at the State Key Laboratory of Cryospheric Science, 
Northwest Institute of Eco-Environment and Resources, CAS. Specifically, we used a laboratory 
standard material (Sigma-Aldrich a-cellulose; Chemical Abstracts Service No. 9004-34-6, 
Germany) to test the effects of measurement conditions on simultaneous measurement of both 
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isotopes. To determine the most suitable reaction temperature for the pyrolysis, we tested the 
results with reaction temperatures ranging from 1000°C to 1440°C. We used a ceramic outer tube, 
a glassy-carbon inner tube, and a graphite electrode during the pyrolysis. We packed 0.16—0.18 
mg Sigma-Aldrich a-cellulose into a silver capsule. High-purity helium (99.9%) was used as the 
carrier gas, at a pressure of 1.2 bar (equivalent to 17 psi). Every sample was run for 800 s. The 
pyrolysis products (CO and H2) were separated using a gas chromatograph under 65°C. 

We defined the optimal temperature of approximately 1400°C for the pyrolysis reaction 
because the temperature produced the smallest deviation from the known isotopic values in the 
laboratory standard material. We heated the samples in a vacuum drying oven at 65°C for 48 h 
before measurements. We then loaded 0.16—0.18 mg a-cellulose of tree rings into a silver capsule. 
Afterwards, the packaged samples were heated again in a vacuum drying oven at 65°C for 24 h 
before the analysis to ensure the samples were dry. For the simultaneous measurement of stable 
isotopes, we set temperatures of 1400°C for the furnace and of 65°C for the gas chromatograph. 
Each tree-ring a-cellulose sample was tested twice. The analytical precision of the standards was 
lower than 0.21%o for 8'8O and 0.10%o for 8'3C. We tested the tree-ring o.-cellulose sample again 
if the standard deviation of the two measurements was higher than 0.3%o for 5!8O and 0.15%o for 
LC. All of the stable isotope results were represented using delta (5) notation and were 
expressed as "%o" values relative to the international standard Vienna Pee Dee Belemnite (VPDB) 
for the carbon isotope and to the Vienna Standard Mean Ocean Water (VSMOW) for the oxygen 
isotope. Hereafter, we refer to these results as the "pyrolysis-1400" series. 


2.2 Correction of mass spectrometer analysis and 61°C corrections 


We calibrated the stable isotope measurements using two methods: the multi-standard calibration 
and the two-point calibration (Evans et al., 2016). The runtime drift correction, amplitude 
correction, and bias correction in mean and deviation were conducted in the two-point calibration 
(Evans et al., 2016). In the multi-standard calibration, we compared the stable isotope values with 
the mean values from four standards with known values. Here, we used the Sigma-Aldrich 
a-cellulose standard and three IAEA (International Atomic Energy Agency) standards 
(IAEA-CH3 standard, benzoic acid IAEA 601, and IAEA-C3 standard) for the calibration (Table 
1). In the two-point calibration, we used two standard cellulose materials (Sigma-Aldrich 
a-cellulose and IAEA-CH3; Table 1), and performed the correction based on the mean difference 
between the measured values and the mean of these two known standard values. 


Table 1 Statics of standard materials during a-cellulose calibration 


SBC 80 
Reference material 
Value (%o) Standard deviation (%o0) Value (%o) Standard deviation (%0) 
IAEA-CH3 —24.72 0.05 32.52° 0.25 
Benzoic acid IAEA 601 —27.50* 0.15 23.30 0.30 
IAEA-C3 —24.91 0.49 32.61° 0.40 
Sigma-Aldrich o-cellulose —24.25* 0.10 28.20° 0.20 


Note: All of the stable isotopes are from the IAEA (International Atomic Energy Agency) products except for the values with a 
superscript letter. All of the work to calibrate the reference materials was conducted by at least two laboratories. è, the value was 
calibrated with USGS24 (83C= —16.05%o (£0.04%o)) and USGS40 (8°C= —26.39%o (£0.09%o)). ”, the value was calibrated from the 
benzoic acid IAEA 601. 


For analysis of long-term trends in the 5!3C, it is necessary to correct for the effect of 
isotopically depleted CO2 in the atmosphere as a consequence of industrialization (Feng and 
Epstein, 1996; McCarroll and Loader, 2004). In addition, tree-ring 5'°C series may incorporate an 
additional trend that reflects a change in the tree's physiological response to increasing 
atmospheric CO2 concentrations (McCarroll et al., 2009; Gagen et al., 2011; Frank et al., 2015). In 
this study, we used the "pin" (pre-industrial) correction method to correct the 5'3C series (8!°Cpin) 
for the effect of the increasing atmospheric CO2 concentrations. The correction was implemented 
in the MATLAB software (McCarroll et al., 2009; McCarroll, 2010). This offered a conservative 
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approach for dealing with the CO2 effect and has been successfully applied in tree-ring ecological 
researches (Andreu-Hayles et al., 2011) and dendroclimatological studies (Hafner et al., 2011; 
McCarroll et al., 2011; Xu et al., 2014; Lavergne et al., 2017). We used the 8'°Cyin to detect the 
climatic signals recorded in tree-ring 8'°C. 


2.3 Climate data 


We used all available climate data (during 1951-2008) from a low-elevation station in Urumqi 
(43°46'48"N, 87°37'12"E; 919 m a.s.l.) located about 70 km northeast of the tree-ring sampling 
site, to explore the climate responses of tree-ring ôC and 8'8O based on combustion and 
pyrolysis methods. These climate records were obtained from the Chinese Meteorological Data 
Center (http://data.cma.cn/), and were also used in our previous study (Xu et al., 2014). The 
homogeneities of these data had been confirmed. The climate variables included mean monthly 
temperature, total monthly precipitation, and mean monthly relative humidity. We also used mean 
monthly values of the standardized precipitation evapotranspiration index (SPED 
(Vicente-Serrano et al., 2010) from the grid point (43°N, 87°E) nearest the sampling site to detect 
the responses of tree-ring 5'°C and 5!80 series to drought. 


2.4 Statistical methods 


We used independent-sample t-test to detect the difference in the mean for 5'3C values determined 
by combustion and pyrolysis methods, and for 8'8O values from the pyrolysis method at 1350°C 
and the simultaneous measurement (pyrolysis at 1400°C). We used one-way analysis of variance 
(ANOVA) to detect differences in 5!°C and 8'8O values between the two methods. We applied 
linear regressions to explore the trends in the C and 6!8O series from 1950 to 2008. To 
determine whether the time series produced by the two methods differ from each other, we 
compared the trends in these series by examining the slope of the regressions. 

We calculated the value of Pearson's correlation coefficient (r) for the relationships between the 
tree-ring isotope data (5'°Cpin and 5!8O) and the climate variables (mean monthly temperature, 
total monthly precipitation, monthly relative humidity, and SPEI) during the common period from 
1951 to 2008. To compare these results with previously published findings (Xu et al., 2014), we 
performed this analysis for the period from September of the previous year (the year before the 
tree-ring formed) to September of the current year. 


3 Results and discussion 


3.1 Dependence of a-cellulose 51°C and 530 on the pyrolysis temperature 


Uncalibrated measured values of tree-ring 5'°C and 8'8O for the Sigma-Aldrich o-cellulose from 
simultaneous measurement both increased with rising furnace temperature, particularly at 
temperatures above 1250°C (Fig. 1). The mean measured 5!8O showed a small difference (2.1%o) 
for the Sigma-Aldrich o-cellulose standard at temperatures from 1000°C to 1250°C. The mean 
measured 5!3C also showed a small difference (0.74%o) among the same temperature ranges. The 
standard deviation of the measurement was relatively large for 5'8O (0.31%0—0.73%o) and the 
number of outliers was relatively high when the furnace temperature was lower than 1250°C (Fig. 
1). These results indicate a relatively low stability of the measurement at temperatures lower than 
1250°C using the TC/EA because the a-cellulose was only partially (about 5/6) converted into 
CO at these temperatures (Leuenberger and Filot, 2007; Young et al., 2011). 

At temperatures from 1300°C to 1440°C, the mean measured 5!8O (without calibration) 
showed a large difference (4.20%o) for the Sigma-Aldrich a-cellulose. The difference between the 
mean measured values was 0.90% for 5!8O and 0.16%o for lC at temperatures of 
1300°C-1350°C (Fig. 1). The difference decreased to 0.49%o for 5!8O but increased to 0.27%o for 
SLC at temperatures of 1400°C-1440°C, which are nearly the acceptable precisions for 5'8O 
(0.30%0) and 5C (0.20%o). The standard deviation was 0.25%o for 5!8O and 0.10%o for 8'3C at 
the temperature of 1400°C (Fig. 1). Thus, high-temperature (1400°C) pyrolysis of tree-ring 
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a-cellulose resulted in a nearly complete conversion of oxygen into CO (Leuenberger and Filot, 
2007; Woodley et al., 2012). Measurement at 1400°C was more stable and precise than 
measurement at lower temperatures, and was able to prolong the lifetime of the furnace compared 
to higher temperatures. In comparing temperature ranges, previous studies have found that the 
reaction temperature of 1400°C appears to be optimal for the simultaneous measurement of 
o.-cellulose 5!8O and 5!3C (e.g., Woodley et al., 2012). 
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Fig. 1 Temperature-dependence of the simultaneously measured (a) 8'8O and (b) 8'°C from the Sigma-Aldrich 
a-cellulose. The box plots represent the measured values, median (horizontal line), minimum and maximum 
values (whiskers), and the 25'* to 75'™ percentiles. The open dots represent outliers (values more than 2 SD 
(standard deviation) from the mean). The measured 8!8O and 8!°C values were not calibrated to the standard 
values (the Vienna Pee Dee Belemnite (VPDB) and the Vienna Standard Mean Ocean Water (VSMOW)) in order 
to exhibit the repeatability (standard errors and outliers) of the measurements. The curve and the equation indicate 
the relationship between the stable isotope value (y) and the pyrolysis temperature (x). 


3.2 Calibration of the simultaneous measurement of °C and 6!8O 


We compared the results using different calibrated methods with data from simultaneous 
measurement (Fig. 2). The 5'°C from simultaneous measurement showed a clear decreasing trend 
in the two-point and multi-standard calibrations (Fig. 2a). Specifically, the multi-standard 
calibrated 5'°C decreased at an annual rate of —0.057%o, whereas the two-point calibrated 85°C 
decreased at an annual rate of —0.041%o. The standard deviation of the multi-standard calibrated 
SLC (1.05%) was higher than that of the two-point calibrated 5!°C (0.81%) for the period 
1950-2008 (Fig. 2a). There were some differences between the lC values from the 
multi-standard calibration and from the combustion method, especially in the late 20" century, 
which resulted from the significantly depleted 8'C in atmospheric CO2 at this time period 
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(Woodley et al., 2012). The multi-standard calibrated 53C from simultaneous measurement had 
an offset (mean offset 0.43%) with the 8'°C from the combustion method, and the offset was 
large at the end of the 20" century (Fig. 2a). This result is similar to the findings from other 
studies (Kndller et al., 2005; Leuenberger and Filot, 2007), which reported an offset between the 
tree-ring 8'°C results measured by the combustion method and the high-temperature pyrolysis 
method (at 1400°C). The multi-standard calibrated °C from simultaneous measurement slightly 
underestimated the low 8'°C values obtained from the combustion method (Figs. 2a and b). These 
results differ from previous findings with "low-temperature pyrolysis" conducted at a temperature 
of 1090°C (Young et al., 2011). Specifically, the multi-standard calibrated 6'°C from 
simultaneous measurement slightly overestimated the low 5!°C values and underestimated the 
high 5!3C values, due to the potential influence of the glassy carbon. 
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Fig. 2 (a) Time-series comparison of !?C values during 1950-2008 from 59 tree-ring o-cellulose samples 
combusted to CO2 in an elemental analyzer Flash EA at 1050°C and trends of 5'C values from two calibrated 
methods (multi-standard calibration and two-point calibration; Evans et al., 2016) for the pyrolysis to CO in the 
TC/EA at 1400°C. The "pin" corrected tree-ring °C values from the combustion method and the two-point 
calibrated pyrolysis method were also exhibited. (b) Slopes and linear regressions for the 5'°C values from the 
combustion method and from the pyrolysis method calibrated by the multi-standard and two-point corrections. (c) 
Time-series comparison of tree-ring 8!8O values from 59 a-cellulose samples from the pyrolysis to CO in the 
TC/EA at 1350°C (pyrolysis-1350) and from the pyrolysis method (simultaneous measurement) at 1400°C 
(pyrolysis-1400). The two-point calibrated correction represents that the 8!8O values from the pyrolysis method 
(simultaneous measurement) at 1400°C (pyrolysis-1400) corrected the offset using the equations from Figure la 
(see the text for details). (d) Slopes and linear regressions for the 5'8O values from the pyrolysis method at 
1350°C and the multi-standard calibrated 8'8O values and two-point calibrated 5'8O values after corrected offset 
from the pyrolysis method (simultaneous measurement) at 1400°C. For figures in the right pane, the black 
diagonal line represents the 1:1 line (y=x) and the dashed line represents the regression line. 


The two-point calibrated 5'°C values from the pyrolysis method and the combustion method 
were very similar. Both 5!°C series showed a decreasing trend with an annual rate of —0.041%o, 
and there was no obvious offset (mean of 0.100%) between them (Fig. 2a). After the "pin" 
correction of these two series, the &!°Cpin values showed no significant trend from 1950 to 2008. 
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The 8'°C value from the combustion method was consistent with the value from the pyrolysis 
method for both the two-point and multi-standard calibrations (R’>0.91; Fig. 2b), and both 
regressions were close to the 1:1 line. Our findings were better than the results reported by 
Woodley et al. (2012), who found large differences between the 8'°C values from the combustion 
and pyrolysis methods after the 20" century. These large differences could have been due to 
potential errors in methods, such as preferential deposition of carbon atoms, exchange with the 
glassy carbon, and formation of carbon residues in the carbon tube because of some uncompleted 
pyrolysis (Young et al., 2011). However, after we adopted the two-point calibrated method with 
two celluloses materials as standards, these errors were eliminated greatly because the two-point 
calibrated 5'°C had a smaller standard deviation and the slope of 8'°C from the two-point 
calibrated method and combustion method was closer to 1 (Fig. 2b). Our results also indicate that 
the two-point calibrated method of Evans et al. (2016) improves the precision of the 8'°C from 
simultaneous measurement (high-temperature pyrolysis) of the two isotopes compared with the 
multi-standard calibrated method. 

There was an offset between the tree-ring 5'8O from the pyrolysis method at 1350°C and from 
simultaneous measurement at 1400°C (Figs. 2c and d). The offset was nearly constant (at about 
—1.52%o) for the two-point calibrated 5'8O during 1950-2008, but varied for the multi-standard 
calibrated 5'8O during some periods (e.g., 1970-1990). The years of 1970-1990 were a relative 
wet period (Xu et al., 2014), and the large variation of 5!8O in this period may result in a large 
offset variation. The 5!8O values from simultaneous measurement were higher than those from the 
pyrolysis method at 1350°C, which may result from the increased reaction temperature (Fig. 1). 
Because of the difference in 5!8O from the pyrolysis method at different temperatures, we applied 
an additional correction to the two-point calibrated 5'8O from simultaneous measurement by 
subtracting the difference (1.67%) between the 5'8O values at 1400°C and 1350°C, which was 
calculated by the equations in Figure la. The offset-corrected 5!8O showed a very similar 
variation to the 5'8O from the pyrolysis method at 1350°C, and both of the two 5!8O series 
showed no significant temporal trends (Figs. 2c and d). In addition, the 5!°O from the pyrolysis 
method at 1350°C and from simultaneous measurement at 1400°C exhibited nearly parallel slopes 
(0.8620 and 0.9125 for the multi-standard and two-point calibrations, respectively), both of which 
were close to the 1:1 line, indicating that simultaneous measurement will not change the 
variability of the 5'8O from the pyrolysis method at 1350°C (Fig. 2d). 


3.3 ANOVA test of the simultaneous and separate measurements 


The mean 8'3C values for the P. schrenkiana samples based on the combustion and two-point 
calibrated simultaneous measurements were very similar (Fig. 3a) and did not differ significantly 
(—22.05%o0 (+0.05%o0) for the combustion method and —22.15%o (+0.10%o0) for simultaneous 
measurement; ANOVA test, P=0.513). The standard deviations of the 59 tree-ring a-cellulose 
samples were also similar (0.83%o and 0.81%, respectively, for the combustion method and 
simultaneous measurement). The range (maximum minus minimum) of the ôC from the 
combustion method was slightly larger than that from simultaneous measurement. The differences 
between the 5°C values from the combustion method and simultaneous measurement were 
normally distributed (Shapiro-Wilk test, P=0.35), with most of the values in the range from 
—0.3%0 to 0.4% (Fig. 3b). 

The mean 5!8O values for the 59 tree-ring o-cellulose samples from the separate pyrolysis 
method at 1350°C and from the temperature-corrected simultaneous measurement at 1400°C were 
also very similar and did not differ significantly (—32.31%0 (+0.25%o) for pyrolysis at 1350°C and 
—32.16%o0 (+0.25%o0) for simultaneous measurement at 1400°C; ANOVA test, P=0.386; Fig. 3c). 
The standard deviations were also similar, with the values of 0.87%o and 0.97%o for pyrolysis at 
1350°C and for simultaneous measurement at 1400°C, respectively. The differences between the 
two methods were normally distributed (Shapiro-Wilk test, P=0.30), and most of the differences 
ranged from —0.5%o to 0.5%o (Fig. 3d). The range of the 5!8O values from simultaneous 
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measurement was larger than that from the pyrolysis method at 1350°C (Fig. 3d). This means that 
after the correction of the mean 8!8O difference caused by the pyrolysis temperature, the means 
and standard deviations of 5'8O from simultaneous isotopic measurement were similar with the 
values obtained from separate isotopic analyses. 
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Fig. 3 (a) Comparison of 8'°C values from 59 tree-ring a-cellulose samples based on combustion and pyrolysis 
methods; (b) the density (width=0.1) distribution and density distribution curve (red line) of the differences 
between paired annual tree-ring 5!°C values determined using the combustion and pyrolysis methods; (c) 
comparison of 8'8O values from 59 tree-ring a-cellulose samples using the pyrolysis method at 1350°C 
(pyrolysis-1350) and the pyrolysis method (simultaneous measurement) at 1400°C (pyrolysis-1400); (d) the 
density (width=0.1) distribution and density distribution curve (red line) of the differences between paired annual 
tree-ring 5'8O values from the pyrolysis method at 1350°C (pyrolysis-1350) and the pyrolysis method 
(simultaneous measurement) at 1400°C (pyrolysis-1400). In the left pane, values are means (purple 
dots)+standard deviations (error bars) for 59 tree-ring a-cellulose samples (yellow dots). Values in the left pane 
did not differ significantly (ANOVA-test) between the two methods. 


The simultaneously measured 8'°C and 5!8O values (pyrolysis method) were typically similar 
to the separately measured values (combustion and pyrolysis methods, respectively), with 
non-significant differences between them. Simultaneous measurement may be considered less 
accurate than separate measurement, especially for the 5'°C determined by the combustion 
method. Young et al. (2011) suggested that if suitable laboratory procedures (such as those 
described in their paper) were adopted, any small loss in accuracy would be outweighed by the 
greater efficiency of simultaneous measurement of two isotopes (i.e., the potential to greatly 
decrease the sample size). For tree rings, the variability between samples from different 
individuals (1%—3%o for 5'°C and 1%o—4%o for 8'8O; Leavitt, 2010) was much greater than the 
analytical uncertainty and the isotopic offset in different measurement methods (—0.3%o0—0.4%o0 for 
LC and —1.0%c-1.5%o for 5'8O; Fig. 3), and even for trees growing in close proximity, the 
difference between individuals can be larger than the difference in the records of climate change 
within an individual tree (McCarroll and Loader, 2004; Leavitt, 2010). Although our research was 
conducted in an arid area, the results are consistent with findings from humid areas (Young et al., 
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2011; Woodley et al., 2012). That is, the 5'°C and 5!8O values from simultaneous measurement 
are similar to those from separate measurement. This indicates that the simultaneously measured 
LC and 8'8O can be applied both in humid and arid areas. In climate reconstructions, the 
precision of the estimated isotopic average values in tree rings depends much more on the number 
of trees in the sampling than on the precision of individual measurements of the isotope values 
(McCarroll and Loader, 2004; Young et al., 2011). Therefore, the simultaneous measurement 
could produce an acceptable regional isotope chronology due to the small decrease of accuracy, 
combined with the similar mean values and trends for tree-ring 5'°C and 5!8O (Figs. 2 and 3). 


3.4 Climate response of tree-ring ðC and 6%O 


The relationships between the stable isotopic values (8'°C and 8!°O) and the four climate 
variables (mean monthly temperature, total monthly precipitation, mean monthly relative 
humidity and SPEI) were very similar, regardless of whether the isotopic values were from 
separate measurement or from simultaneous measurement (Fig. 4). The climate response inferred 
from the combustion 6'Cpin was significantly (P<0.01) positively correlated with mean 
temperature in July (r=0.33), and negatively correlated with relative humidity in July, SPEI in 
July, and SPEI of the whole growing season from May to August(r= —0.38, P<0.01; Fig. 4). The 
climate response results for the 5!°Cpin from the pyrolysis method suggest that 5'C was 
significantly positively correlated with the overall mean temperature of the whole growing season 
(r=0.40, P<0.01), and was significantly (P<0.01) negatively correlated with precipitation and 
SPEI during the whole growing season (Fig. 4). The &!°Cpin from the pyrolysis method did not 
alter the pattern of 5'°C-climate response and kept the climatic signals. Furthermore, it reveals 
that some significant climatic signals were not detected by the combustion method (e.g., negative 
and significant (P<0.01) correlations for temperature in December of the previous year, and for 
SPEI in May of the current year; Fig. 4). The 8'C,in from the pyrolysis method enhanced the 
climatic signal strength of overall mean temperature, precipitation, relative humidity, and SPEI 
for the whole growing season (Fig. 4). These results indicate that the 5'°Cyin from the pyrolysis 
method does not greatly alter the interpretation of the climate response of tree-ring 61°C. 

Tree-ring 8!8O values from the pyrolysis method at 1350°C and from simultaneous 
measurement at 1400°C were significantly (P<0.01) positively correlated with mean monthly 
temperature from May to August (7>0.36) and the overall average temperature of the whole 
growing season (r=0.54). However, they were significantly (P<0.01) negatively correlated with 
precipitation in July, relative humidity from May to July, and SPEI in May and July (Fig. 4). Both 
of the two tree-ring 5'8O series (from the pyrolysis method at 1350°C and from simultaneous 
measurement at 1400°C) showed significant correlations (P<0.01) with the overall mean 
temperature (7=0.54 and r=0.53, respectively), precipitation (r= —0.42 and r= —0.44, respectively), 
relative humidity (r= —0.49 and r= —0.49, respectively) and SPEI (r= —0.51 and r= —0.54, 
respectively) for the whole growing season (Fig. 4). The tree-ring 5'°O from simultaneous 
measurement, therefore, did not alter the dominant climatic signals in the tree-ring 5!8O series, 
and clearly reflected both temperature and moisture conditions during the growing season (Xu et 
al., 2014). As was the case for 5'°Cpin, the 8180 from simultaneous measurement reveals that some 
significant (P<0.01) climatic signals were not detected in separate measurement (e.g., 
precipitation in April of the current year; Fig. 4). 

Because the climate response analyses or climate reconstruction would re-scale the series of 
tree-ring stable isotopes, the climatic signals from the pyrolysis and combustion methods would 
change little (Young et al., 2011). Taken together, our results indicate that the climate response of 
tree-ring 5°C and 8'8O obtained from simultaneous measurement is the same as that from 
separate measurement, and therefore simultaneous measurement can be used as an alternative 
method for detecting the sable isotopic values (C and 6!8O). These good results can be 
achieved due to the low loss of accuracy and the similar regional mean values and trends provided 
by the simultaneous measurement (Figs. 2 and 3). Thus, for analysis of a regional tree-ring 
isotope network, the simultaneous measurement would greatly improve the number of sampling 
points that can be analyzed, and will therefore improve the regional coverage. 
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Fig. 4 Climate responses of tree-ring 8'Cpin and 8!8O series determined from separate and simultaneous 
measurements. For simultaneous measurement, the stable isotope values were calibrated with the two-point 
method (Evans et al., 2016). The tree-ring &!?Cpin was the 5!5C corrected by using the method of "pin" correction 
(McCarroll et al., 2009). The size and color of the circles indicate the strength and significance of the relationship 
(Pearson's correlation coefficient), respectively. The red and blue circles indicate the Pearson's correlations are 
statistically significant (P<0.01). Months with a "p" prefix represent months in the previous year. MJJA represents 
the mean value from May to August in the growing season of the current year. Temp, mean monthly temperature; 
Pre, total monthly precipitation; RH, mean monthly relative humidity; SPEI, monthly mean standardized 
precipitation-evapotranspiration index. 


4 Conclusions 


In this study, we compared the 8'°C and 5!8O values obtained from separate and simultaneous 
measurements to determine the suitability of simultaneous measurement using the 
high-temperature pyrolysis method in an arid region of northwest China. We found that the 6!°C 
values from the combustion method were very similar to those from simultaneous measurement 
(high-temperature pyrolysis method). The two-point calibration method improved the results for 
both 53°C and §!8O, and the mean value, standard deviation, and trend of the calibrated §'C 
series from the pyrolysis method (simultaneous measurement) were similar to those obtained 
from the combustion method (separate measurement). The simultaneously measured 5!8O values 
from the pyrolysis method at 1400°C showed a nearly constant offset compared with the values 
from the pyrolysis method at 1350°C because 5!8O depends on the reaction temperature. However, 
the variability of the 5'°O from simultaneous measurement was similar to that from separate 
measurement in the time series. The correlation of the stable isotope values with temperature and 
moisture variables indicates that simultaneous measurement will not change the inferred climatic 
signals. Based on the results mentioned above, we conclude that simultaneous measurement of the 
8C and 5!80 values in tree-ring a-cellulose samples through high-temperature pyrolysis method 
would produce acceptable and reliable stable isotope series in both arid and humid areas. The 
simultaneous isotopic measurement can greatly reduce the cost and time requirement compared 
with separate isotopic measurement. 


XU Guobao et al.: Application and verification of simultaneous determination of cellulose 8!3C and... 


Acknowledgments 


This research was funded by the National Natural Science Foundation of China (41501049, 41571196), the 
Self-determination Project of the State Key Laboratory of Cryospheric Sciences (SKLCS-ZZ-2018), the "Light of 
West China" Program of the Chinese Academy of Sciences, the Youth Innovation Promotion Association, Chinese 
Academy of Sciences (2016372), the Chinese Scholarship Council (201704910171), and the Fundamental 
Research Funds for the Central Universities (GK201801007). The authors thank Prof. Michael EVANS from the 
University of Maryland, USA, and Prof. Iain ROBERTSON from Swansea University, UK for providing the 
calibration code and advice on calibration materials. We also thank Prof. LIU Weiguo from Institute of Earth 
Environment, Chinese Academy of Sciences, Miss ZHANG Pingyu from Lanzhou University, and Miss ZHANG 
Li from Institute of Botany, Chinese Academy of Sciences, for calibrating the standard materials. Special thanks 
to the anonymous reviewers and editors for their constructive comments on this manuscript. 


References 


Altieri S, Mereu S, Cherubini P, et al. 2015. Tree-ring carbon and oxygen isotopes indicate different water use strategies in three 
mediterranean shrubs at Capo Caccia (Sardinia, Italy). Trees-Structure and Function, 29(5): 1593-1603. 

Andreu-Hayles L, Planells O, Gutiérrez E, et al. 2011. Long tree-ring chronologies reveal 20" century increases in water-use 
efficiency but no enhancement of tree growth at five Iberian pine forests. Global Change Biology, 17(6): 2095-2112. 

Babst F, Alexander M R, Szejner P, et al. 2014. A tree-ring perspective on the terrestrial carbon cycle. Oecologia, 176(2): 
307-322. 

Dawson T E, Siegwolf R T W. 2007. Using stable isotopes as indicators, tracers, and recorders of ecological change: Some 
context and background. Terrestrial Ecology, 1: 1, 3-18. 

Evans M N, Selmer K J, Breeden B T, et al. 2016. Correction algorithm for online continuous flow 8'°C and 5!8O carbonate and 
cellulose stable isotope analyses. Geochemistry, Geophysics, Geosystems, 17(9): 3580-3588. 

Farquhar G D, O'Leary M H, Berry J A. 1982. On the relationship between carbon isotope discrimination and the intercellular 
carbon dioxide concentration in leaves. Australian Journal of Plant Physiology, 9: 121—137. 

Feng X, Epstein S. 1996. Climatic trends from isotopic records of tree rings: The past 100-200 years. Climatic Change, 33(4): 
551-562. 

Frank D C, Poulter B, Saurer M, et al. 2015. Water-use efficiency and transpiration across European forests during the 
anthropocene. Nature Climate Change, 5(6): 579-583. 

Gagen M, Finsinger W, Wagner-Cremer F, et al. 2011. Evidence of changing intrinsic water-use efficiency under rising 
atmospheric CO concentrations in boreal fennoscandia from subfossil leaves and tree ring 8'°C ratios. Global Change 
Biology, 17(2): 1064-1072. 

Hafner P, Robertson I, McCarroll D, et al. 2011. Climate signals in the ring widths and stable carbon, hydrogen and oxygen 
isotopic composition of Larix decidua growing at the forest limit in the southeastern European Alps. Trees-Structure and 
Function, 25(6): 1141-1154. 

Keel S G Joos F, Spahni R, et al. 2016. Simulating oxygen isotope ratios in tree ring cellulose using a dynamic global 
vegetation model. Biogeosciences, 13(13): 3869-3886. 

Keller K M, Lienert S, Bozbiyik A, et al. 2017. 20" century changes in carbon isotopes and water-use efficiency: 
tree-ring-based evaluation of the CLM4.5 and LPX-Bern models. Biogeosciences, 14(10): 2641-2673. 

Kłusek M, Pawelczyk S. 2014. Stable carbon isotope analysis of subfossil wood from Austrian Alps. Geochronometria, 41(4): 
400-408. 

Knöller K, Boettger T, Weise S M, et al. 2005. Carbon isotope analyses of cellulose using two different on-line techniques 
(elemental analysis and high-temperature pyrolysis)—a comparison. Rapid Communications in Mass Spectrometry, 19(3): 
343-348. 

Laumer W, Andreu L, Helle G et al. 2009. A novel approach for the homogenization of cellulose to use micro-amounts for 
stable isotope analyses. Rapid Communications in Mass Spectrometry, 23(13): 1934—1940. 

Lavergne A, Daux V, Villalba R, et al. 2017. Improvement of isotope-based climate reconstructions in patagonia through a 
better understanding of climate influences on isotopic fractionation in tree rings. Earth and Planetary Science Letters, 459: 
372-380. 

Leavitt S W. 2010. Tree-ring C-H-O isotope variability and sampling. Science of the Total Environment, 408(22): 5244-5253. 

Leavitt S W, Treydte K, Liu Y. 2010. Environment in time and space: opportunities from tree-ring isotope networks. In: West J, 
Bowen G, Dawson T, et al. Isoscapes. Dordrecht: Springer, 113—135. 

Leuenberger M C, Filot M S. 2007. Temperature dependencies of high-temperature reduction on conversion products and their 


isotopic signatures. Rapid Communications in Mass Spectrometry, 21(10): 1587—1598. 


JOURNAL OF ARID LAND 


Liu X, An W, Leavitt S W, et al. 2014. Recent strengthening of correlations between tree-ring 5'°C and 5'%O in mesic western 
China: Implications to climatic reconstruction and physiological responses. Global and Planetary Change, 113: 23-33. 

Liu Y, Cobb K M, Song H, et al. 2017. Recent enhancement of central pacific El niño variability relative to last eight centuries. 
Nature Communications, 8: 15386, doi: 10.1038/ncomms15386. 

Loader N J, Robertson I, Switsur V R, et al. 1997. An improved technique for the batch processing of small wholewood samples 
to a-cellulose. Chemical Geology, 136(3—4): 313-317. 

Loader N J, Street-Perrott F A, Daley T J, et al. 2015. Simultaneous determination of stable carbon, oxygen, and hydrogen 
isotopes in cellulose. Analytical Chemistry, 87(1): 376-380. 

McCarroll D, Loader N J. 2004. Stable isotopes in tree rings. Quaternary Science Reviews, 23(7—8): 771-801. 

McCarroll D, Gagen M H, Loader N J, et al. 2009. Correction of tree ring stable carbon isotope chronologies for changes in the 
carbon dioxide content of the atmosphere. Geochimica et Cosmochimica Acta, 73(6): 1539-1547. 

McCarroll D. 2010. Erratum to D. McCarroll, M. Gagen, N. Loader, I. Robertson, K. Anchukaitis, S. Los, G Young, R. 
Jalkanen, A. Kirchhefer, J. Waterhouse (2009), "Correction of tree ring stable carbon isotope chronologies for changes in the 
carbon dioxide content of the atmosphere", Geochimica et Cosmochimica Acta 73, 1539-1547. Geochimica et Cosmochimica 
Acta, 74(10): 3040, doi: 10.1016/j.gca.2009.12.031. 

McCarroll D, Tuovinen M, Campbell R, et al. 2011. A critical evaluation of multi-proxy dendroclimatology in northern Finland. 
Journal of Quaternary Science, 26(1): 7—14. 

Roden J S, Lin G, Ehleringer J R. 2000. A mechanistic model for interpretation of hydrogen and oxygen isotope ratios in 
tree-ring cellulose. Geochimica et Cosmochimica Acta, 64(1): 21-35. 

Roden J S, Farquhar G D. 2012. A controlled test of the dual-isotope approach for the interpretation of stable carbon and oxygen 
isotope ratio variation in tree rings. Tree Physiology, 32(4): 490-503. 

Saurer M, Siegenthaler U, Schweingruber F. 1995. The climate-carbon isotope relationship in tree rings and the significance of 
site conditions. Tellus B, 47(3): 320-330. 

Saurer M, Spahni R, Frank D C, et al. 2014. Spatial variability and temporal trends in water-use efficiency of European forests. 
Global Change Biology, 20(12): 3700-3712. 

Scheidegger Y, Saurer M, Bahn M, et al. 2000. Linking stable oxygen and carbon isotopes with stomatal conductance and 
photosynthetic capacity: a conceptual model. Oecologia, 125(3): 350-357. 

Schleser G H, Helle G, Lücke A, et al. 1999. Isotope signals as climate proxies: the role of transfer functions in the study of 
terrestrial archives. Quaternary Science Reviews, 18(7): 927—943. 

Stokes M A, Smiley T L. 1968. An Introduction to Tree-ring Dating. Chicago: University of Chicago Press, 1-61. 

Treydte K, Frank D, Esper J, et al. 2007. Signal strength and climate calibration of a European tree-ring isotope network. 
Geophysical Research Letters, 34(24): L24302, doi: 10.1029/2007 1031106. 

Vicente-Serrano S M, Begueria S, Lépez-Moreno J I. 2010. A multiscalar drought index sensitive to global warming: the 
standardized precipitation evapotranspiration index. Journal of Climate, 23(7): 1696-1718. 

Wang W, Liu X, Xu G, et al. 2013. Moisture variations over the past millennium characterized by qaidam basin tree-ring 5'*O. 
Chinese Science Bulletin, 58(32): 3956-3961. 

Werner C, Schnyder H, Cuntz M, et al. 2012. Progress and challenges in using stable isotopes to trace plant carbon and water 
relations across scales. Biogeosciences, 9(8): 3083-3111. 

Woodley E J, Loader N J, McCarroll D, et al. 2012. High-temperature pyrolysis/gas chromatography/isotope ratio mass 
spectrometry: simultaneous measurement of the stable isotopes of oxygen and carbon in cellulose. Rapid Communications in 
Mass Spectrometry, 26(2): 109-114. 

Xu G. 2014. Climatic significance of stable oxygen (5'*O) in tree-ring in north part of Xinjiang Uygur Autonomous Region. 
PhD Dissertation. Beijing: University of Chinese Academy of Sciences. (in Chinese) 

Xu G, Liu X, Qin D, et al. 2014. Drought history inferred from tree ring '°C and 65'8O in the central Tianshan Mountains of 
China and linkage with the North Atlantic Oscillation. Theoretical and Applied Climatology, 116(3—4): 385—401. 

Young G H F, Loader N J, McCarroll D. 2011. A large scale comparative study of stable carbon isotope ratios determined using 
on-line combustion and low-temperature pyrolysis techniques. Palaeogeography, Palaeoclimatology, Palaeoecology, 300(1-4): 
23-28. 


